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Abstract 

We study multi-soliton solutions of the four-dimensional SU(N) Skyrme model by combining the 
hedgehog ansatz for SU(N) based on the harmonic maps of into CP^~^ and a geometrical trick 
which allows to analyze explicitly finite-volume effects without breaking the relevant symmetries of 
the ansatz. The geometric set-up allows to introduce a parameter which is related to the’t Hooft 
coupling of a suitable large N limit, in which TV —^ 00 and the curvature of the background metric 
approaches zero, in such a way that their product is constant. The relevance of such a parameter 
to the physics of the system is pointed out. In particular, we discuss how the discrete symmetries 
of the configurations depend on it. 


1 Introduction 

One of the most intriguing theoretical results in Quantum Field Theory (QFT henceforth) has been 
the realization that fermions can emerge out of purely bosonic Lagrangians as solitonic excitations 
(for a detailed review see [T]). The clearest demonstration that the importance of this result goes far 
beyond pure theoretical physics is given by the Skyrme theory [2] which is one of the most important 
models of nuclear and particle physics. The Skyrme term [2] allows the existence of static soliton 
solutions with finite energy, called Skyrmions (see laiais]) describing fermionic degrees of freedom 
(see [61 [3 El [g ini El 12] and references therein). The wide range of applications of the theory in 
other areas (such as astrophysics, Bose-Einstein condensates, nematic liquids, multi-ferric materials, 
chiral magnets and condensed matter physics in general [laEiiisiiiiiiiziiisiiisiiiniiziiiia) is well 
recognized by now. It is also worth to emphasize that the Skyrme model appears in a very natural 
way in the context of the AdS/CFT correspondence |23j . 

From the point of view of nuclear physics, it is very important to have analytic tools allowing to 
analyse the Skyrme model when interacting Skyrmions are present within bounded regions. This case 
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is relevant whenever one wants to take into account the effects of the finite size on the topological 
properties of the Skyrmions themselves. It is commonly believed that many-nucleons systems (such 
as the ones occurring in nuclear pasta |24l [25] ; for a review see [26] ) are completely out of reach of the 
analytical techniques provided by soliton theory already in the SU{2) case (while, at a first glance, 
the SU{N) case is even worse). In particular, the task to compute physical parameters of such multi¬ 
nucleon systems with analytic multi-Skyrmionic configurations is believed to be completely hopeless. 

Recently, the generalized hedgehog ansatz in the SU{2) case, introduced in [271[28l[29|, allowed 
the construction of the first multi-Skyrmions at finite volume: namely, the first exact solutions of the 
Skyrme model representing interacting elementary Skyrmions with a non-trivial winding number, in 
which finite-volume effects can be explicit taken into account (arriving at a good prediction for the 
compression modulus), was obtained |30)^l . The way to do this is to write the system in a modified 
“cylinder-like” metric whose curvature is parametrized by a length Rq. Then multi-Skyrmionic con¬ 
figurations look like necklaces of elementary Skyrmions interacting in bounded tube-shaped regions. 
The ground state of such multi-Skyrmions has the remarkable property that, although the BPS bound 
in terms of the winding cannot be saturated, a new topological charge exists which leads to a different 
BPS bound, which can instead be saturated. 

In this paper, we consider the SU{N) case. A very powerful technique to construct multi- 
Skyrmionic configurations in unbounded regions is given by an ansatz for SU(N) Skyrmions introduced 
in m and |34] . based on harmonic maps of into CP^ which allowed the construction of many 
interesting numerical multi-Skyrmionic configurations. Here we will exploit the fact that the metric 
we use is spherically symmetric, just like flat space, so the harmonic map ansatz can be used also in 
the present case without modifications. We shall see that in the new metric the equations simplify 
with respect to those studied in [33] and [M]; in particular they become autonomous, thus techniques 
from dynamical systems theory become available. Thanks to the choice of the background geometry, 
a novel type of large N limit becomes possible in which N is large and the curvature is small in such 
a way that their product is constant. In this way one can see that both the effects of the curvature 
become negligible and the field equations remain autonomous (so that, in such a limit, one can kill 
the curvature of the metric but keeping all the advantages of the technique of [30] 1. It is also worth 
to emphasize that it is precisely the large N limit which discloses in the clearest possible way the role 
of the Skyrme model as low energy limit of QCD. A non constructive proof of the existence of such 
nontrivial solutions is also provided and a few numerical solutions are exhibited as well. 

This paper is organized as follows: in the second section, the SU{2) case treated in [30] will be 
briefly reviewed in order to set the stage. In the third section, the ansatz for the Skyrmions will be 
described. In the fourth section, the general equations of motion are written down and the flat-large 
N limit is discussed. In the fifth section, we study the system using techniques from the theory 
of dynamical systems; in particular, the stability of the fixed points in the 517(3) case is analyzed, 
and their dependence on the geometric ’t Hooft parameter is discussed. In the sixth section, some 
nontrivial numerical solution of the equations of motion in the 517(4) case will be studied. In the 

^Using similar techniques (see m and [32]), intriguing cosmological properties of the Skyrme model have been 
disclosed. 
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seventh section, some conclusions will be drawn. In the appendix, it is proved on general grounds 
that the equations of motion do have nontrivial solutions which can be interpreted as genuine SU{N) 
Skyrmions. 


2 Generalized hedgehog at finite volume: the SU{2) case 

In this section, the construction of analytic multi-Skyrmionic configurations in the SU(2) case [30] 
will be shortly reviewed. The action of the four-dimensional Skyrme model is given by 

Ssk = j I d^Xy/^Tr > ( 1 ) 

K >0 , A > 0 , Rf,:= U-^Vf,U = := [R^, R^] , 


where the Planck constant and the speed of light have been set to 1, and K and A are the coupling 
constants. The P’s are the generators of the flavor group, which in this section is SU{2). Notice that 
we allowed for curved metrics. The coupling constants K and A are related to the couplings F.,^ and 
e used in [oj^ by 


K = -F‘^ 

4 ’f’ 


K\ = ^ 


The non-linear sigma model term of the Skyrme action is necessary to take into account pions. The 
second term is the only covariant term leading to a well-defined Hamiltonian formalism in time which 
supports the existence of Skyrmions. 

The field equations following from the above action are 


+ = 0 . 

The following standard parametrization of the S'C/(2)-valued scalar U{x^) will be adopted 


( 2 ) 


= y°(x^)i2 ± Y\x>^)ti , (yo) + Y% = 1, 


(3) 


where I 2 is the 2 x 2 identity matrix; to describe a spherically symmetric field configuration we use 
the hedgehog ansatz |28l [30| : 


Y^ = cos a , Y'’ = n^ sin a , a = a{x, t) , (4) 

= sin 9 cos (/? , = sin 9 sin 99 , = cos 9 . (5) 

In order to mimic finite-volume effects without loosing the nice properties of the hedgehog ansatz 
^Experimentally, IT = 186 MeV , e = 5.45. 
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we will consider the following curved 


ds'^ = —dt^ + dx'^ + R^idO'^ + (sin0)^d(/?^) , — ^ < x < ^ , (6) 

0 < 0 < TT , 0 < (/j < 27r , (7) 

where L is the length of the x—interval. The total volume of space is y = AttRqL. This geometry 
describes three-dimensional cylinders whose sections are 5^ spheres, so that parameter Rq plays the 
role of the (finite) diameter of the transverse sections of the tube. The fact that this parameter 
replaces the radial variable r in the metric also leads, as we will see, to considerable simplification 
of the equations of motion, even allowing to find exact solutions in the SU{2) case. Moreover, the 
curvature of this metric is proportional to l/i?Q. As it will be explained in the next sections, the 
explicit presence of this parameter in Eq. ([6|) together with the N of SU{N) allows to define a 
smooth flat limit in which Rq ^ oo and so all the effects of the curvature disappear (however, the 
global topology of space remains cylindrical even in the flat limit and so it differs from the trivial 
topology of flat static unbounded Skyrmions). Thus, in a sense, the above metric is introduced just as 
a “regulator” whose local effects can be removed at the end. In the present context, “flat limit” really 
means 

> 1 , ( 8 ) 

so that, from the practical point of view, already when Rq is around 100 fm all the effects of the 
curvature are negligible and, consequently, even in the flat limit in Eq. ([8]) finite volume effects will 
not disappear. It is also worth to emphasize that the well known result that elementary Skyrmions 
should be quantized as Fermions (which originally was derived on flat spaces) has been extended to 
space-times with compact orientable three-dimensional spatial sections in [38] (and the metric in Eq. 
(|6|) belongs to this class). 

Since, at the end, we will be interested in the flat limit in Eq. Q, one may wonder whether it 
would be possible to start from the very beginning with a flat metric. In fact, as it has been shown 
in m, the background metric in Eq. (|6|) is a very suitable tool to take into account finite volume 
effects (since the total spatial volume is finite) without breaking relevant symmetries of the hedgehog 
ansatz, with the additional advantage of simplifying the field equations. A further relevant advantage 
of the above background metric is that, unlike what happens in the usual unbounded case, it allows 
to define in a very transparent way a smooth large N limit of the SU{N) Skyrmions. Therefore, it is 
much more convenient to analyze the Skyrme theory first within the background metric in Eq. (jUj). 
and take the flat limit only later. 

The effectiveness of such a choice for the metric is also shown by the results in [3S] and [1D| in 
which it has been shown that, unlike what happens in flat space, the equations for the Yang-Mills- 

®The simplest choice of a bounded spatial metric is S®, but this has already been considered in | 35[I36| . The geometry 
in Eq. was considered in but with a motivation different from the analysis of finite-volume effects. Consequently, 
the main results obtained in |30| (the derivation of both a novel BPS bound which can be saturated and analytic multi- 
Skyrmions at finite volume and an explicit formula for the compression modulus in good agreement with experiments) 
are indeed novel. 
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Higgs system (in the sector with non-vanishing non-Abelian electric and magnetic charges) possess 
analytic solutions even in the case in which the Higgs coupling is non-zero. 

With the above ansatz the Skyrme field equations reduce in the static case to the following scalar 
differential equation for the Skyrmion profile a |28| : 


. 2A . 2 , 

1 -|- sm a ] a — 




Cq / .ICQ 

The winding number W for such a configuration reads: 


■ 2 sin a 
-^0 


= 0 


(9) 


=-/ [a sin"^ a ] dx 

TT 


W = j {U-^d^U) {U-^djU) {U-^dkU) 

In the present case, the natural boundary conditions correspond to the choice: 


( 10 ) 


a , nEZ. (11) 


and with these boundary conditions the winding number takes the integer value n. These boundary 
conditions are unique in that they ensure U{—^) = (—1)”[/(-|), which correspond to bosonic and 
fermionic states for even and odd n, respectively. 

Smooth solutions exist for any n satisfying the above boundary conditions for a finite range 
(—L/2,L/2). In particular multi-soliton solutions exist, which represent Skyrmions with winding 
number n living in a finite spatial volume V = AttRqL. It is worth to remark that the large n limit 
in the present context is quite natural in order to consider thermodynamical properties of the multi- 
Skyrmions system and n is the baryon number: obviously, a thermodynamical analysis only makes 
sense in the cases in which the number of particles is very large. 


3 The hedgehog ansatz for SU(N) 

Now we switch to the SU{N) case with generic N. Our analysis will be based on the techniques 
introduced in |33ll34j . which we shall now briefly describe. The hedgehog ansatz for SU{N) spherically 
symmetric Skyrmions living in flat Minkowski metric 

= —dt^ + dx^ + x^{d9‘^ + (sin 9)'^dip‘^) , (12) 

is based on a suitable family of projectors from 3“^ into CP^~^ (see m ). Such projectors P{V) can 
be written as 


P{V) 


|H|2 


(13) 


where |H) is an A^-component complex vector of two complex variables ^ and ^ which locally parametrize 
S'^, where ^ = tan(0/2) and \V\^ = iy\V) . The first P is constructed with a holomorphic 
|H) = |/(^)) while the Gram-Schmidt procedure gives rise to the others. Indeed, P+ can be found by 
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its action on any vector |r;) G [H] as 


P+\v) =d^\v) -lu) (14) 

\v\^ 

Consequently, the other vectors P^\v) are determined inductively: P^\v) = P^{P^~^\v) ). 

Thus the operator P^ corresponding to the family of vectors \V) = |14) = 1-? + /) (for / = f{Cj) 
reads 

Pfc = A: = 0,...,1V-1, (15) 

where, due to the orthogonality of the projectors, we have Ylk=o 

Due to the holomorphy of /, the following identities of the above dehned vectors can be proved 

m- 


{Vk\Vi)=0, 
di\Vk) = 


k^l, 


Af-ll 


d, 


\Vk-i 

\Vk-i 


\Vk) 

\Vk-i\ 


(16) 

(17) 


For SU{N) the components of iVjv-i) = 7+” |/), up to an irrelevant overall factor which cancels in 
the projector, are functions of only 

The SU {N) hedgehog ansatz defined in [33l [M] reads then 


U = exp < iaQ ( Pq — 


I 

N 


+ iai Pi - 


N 


H-+ iaN-2 1 Pn-2 — 


N 


= + AiiPii) e-*“fo^(/ + ^iPi) ... e-“^- 2 /^(/ + AN-2PN-2) 


(18) 


where we defined A^ = e*"*" — 1. Such an ansatz involves the introduction of — 1 projectors and of 
— 1 profile functions = ak{x), k = 0,..., N — 2] . Note that the projector P/v-i is not included 
in the above formula since it is a linear combination of the others. 

One of the main results in [33] and |33| has been to show that the above ansatz in Eq. (jlSp . 
when inserted into the full Skyrme field equations in the case in which the background metric is the 
standard flat metric in spherical coordinates in Eq. (d, gives rise to a consistent set of N" — 1 coupled 
non-linear differential equations for the A^ — 1 profiles a*. Moreover, such field equations can also be 
derived as stationary equations for the energy functional with respect to variations of the profiles. A 
close inspection of the computations in |34j shows that the main requirement in order for the ansatz 
to work is the 50(3) invariance of the background metric. This then suggests that the above ansatz 
may also work in the finite-volume metric in Eq. (|6|) adopted in [30|. This is what we show in the 
next section. 
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4 Hedgehog ansatz for SU(N) at finite volume 

In this section we switch back to the finite-volume <S'0(3)-invariant metric in Eq. ([B]). A direct 
computation shows that the ansatz in Eq. (|18p . when inserted into the full Skyrme field equations 
in this metric gives rise to a consistent system of iV — 1 coupled autonomous non-linear differential 
equations for the N — 1 profiles a* which is simpler than the flat non-autonomous system analyzed 
in [34] as it will be now discussed. Eurthermore, also in the present case the field equations can be 
derived as stationary equations for the energy functional with respect to variations of the a^’s. It is 
worth emphasizing here that the explicit presence of the parameters N and Rq^L allows to consider 
a flat large N limit in which the curvature is negligible and the product NR^'^L remains constant. 
Hence, the present formalism is also relevant for people only interested in the Skyrme model on flat 
space-times. It is important to notice that in such a flat limit, while the effects of the curvature of the 
metric in Eq.® disappear, the effects of the cylindrical topology do not. 

From now on we shall use the following energy [E] and length [L] units: 

[E] = ^^l GeV , (19) 

[L] = ^ « 0.6 fm ^ (20) 

K = 2 , A = 1 . (21) 

Following the same steps of [34] in the new metric, one arrives at the following expressions for the 
winding number, the total energy of the SU{N) hedgehog and the field equations respectively: 

N-2 

i=0 

\2 

Oik-1 ) Dk 

(23) 
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0 — 7T— ^ (1 + i) Fi + ^ ^ Fi H ^^2-(1 “ cos {Fi)) Fi 

k=0 i=k 
2 


N 


+ 


i=0 

{l + l){N-l-l) 


Rl 


2Rl 


sm(F,)(0) + ' ^^ sm(F) 


iq 


+ 


{l + lf{N -l-lf 

K 

(/ + l)(Ar-/-i)sm(F;) 

2Rt 


(1 - cos (Fi))sin(F/) 


[I [N - 0 (1 - cos (Fz_i)) + {1 + 2){N-l-2){l- cos (F;+i))] , 

(24) 


where we introduced the quantities 


Fk = otk — Oik+i , k = 0,... ,N — 3, Fn_2 = oin-2, 


(25) 


Fi = dFi/dx, and we defined 


Dk = k {N - k) {I - cos {ok - Ok-i)) ■ 


(26) 


Notice that the winding number gets two different contributions (while in the Minkowski case it 
just gets one), as a consequence of the different topology of the cylinder-like metric we are using. 
Another important point is that when the parameter Rq is very large the field equations (I24p tend 
to linear equations, on the other hand when Rq is small the nonlinear terms become very important. 
Thus, this parameter (more precisely, the quantity l/i?o) controls the nonlinearity of the theory, and 
thus it plays the role of an additional coupling constant. All the analyses we shall perform in the 
following will confirm the fact that the physics of the model crucially depends on it. More precisely, 
both the numerical analysis of Sect. [6] and the existence theorem of the appendix tell us that the 
global properties of the solutions of the system depend on the dimensionless parameter Rq/L, which 
controls the shape of the cylinder. 

By comparing the present total energy in Eq. ()23p and the field equations in Eq. (I24p with the 
corresponding expressions in [31], the first of the advantages mentioned above in working within the 
metric in Eq. ([6|) is apparent. Namely, the field equations become an autonomous system which can, 
as such, be analyzed with the powerful tools of dynamical systems theory. 

In the special case in which all the profiles are equal, i.e. Fi = F, V Z all the equations of the above 
system become proportional and equivalent to: 


F 



(1 — cosF) \ 

) 


sinE 



2(1 — cosE)\ 

M y. 


= 0 


(27) 


Upon setting F = 2q;, this is precisely Eq. ([9|), valid in the 517(2) case. This is to be expected since 
setting all the profiles to be equal corresponds to looking for solutions which are embeddings of the 
517(2) ones in SU{N). 

In the next sections we shall see numerically in some cases that the above system of equations 
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does admit nontrivial solutions, i.e. solutions whose profiles are not equal nor proportional, while in 
the appendix we shall prove analytically the existence of such solutions. Such solutions have winding 
numbers which in general are nonvanishing, therefore can be interpreted as genuine SU{N) multi- 
Skyrmion configurations. 


4.1 The flat-large N limit 

As said, the formalism developed so far strongly relies on the usage of a spacetime whose spatial 
sections have a finite radius Rq, and whose curvature goes to zero as iio cc. In this section we will 
show that there is a natural way to accommodate this flat limit, which requires to take the large N 
limit as well. The latter is a common tool in Quantum Field Theory, see e.g. |42] for a review. 

The field equations in Eq. (I24p can be rewritten as follows 


N-2 


I N-2 


0 = - ( -2(1 + 1) (1 + i) F, + 2iV F, I + 


i=0 


k=0 i=k 


{N - 1) F2 


(1 - cos (F/)) Fi 


+ sm 


{N - 1) 2F2 

{N -I-if {I+ iy 

{N -if 


iFi)( 


^,1 




\ / (^ - 1) 
(1 - cos (Fi)) sin (Fi) 


+ 


{N-l-l) (/ + l)sin(Fz) 


(iV-1) 


2R^ 


(1 - cos (F,_i)) + (/ + 2) ^ (1 - cos {Fl+^)) 


(N-l) 


(N-l) 


(28) 


where we introduced the effective radius 


p2 

R^= ^0 


iV- 1 


(29) 


The nonlinear part of the field equations (|28p behaves smoothly in the large N limit provided also 
Ro is large, so that R is kept constant. More precisely: 

r2 

N ^ oo , Fn —^ oo I lim -— = R? = finite . (30) 

Hence, the proper way to consider the large N limit is to simultaneously consider the flat limit Fq —?> oo 
in such a way that the parameter F^ in Eq. (I29p remains finite. Since in the large N limit in Eq. 
(|30l) the nonlinear part of the field equations Eq. (1281) does not depend on N and Fq separately but 
only on the effective radius F defined in Eq. (|2^ . the quantity 1/F^ = {N — 1)/Fq plays the role 
of a geometric (it defines a length scale) ’t Hooft coupling. Physically, the parameter R^ in Eq. (j29p 
represents the “effective area” available for each Skyrmion within each section of the tube. 

In conclusion, one can see the metric defined in Eq. ([6|) just as a technical device in order to analyze 
multi-Skyrmionic configurations since, at the end, one can turn off all the effects of the curvature 
(keeping the cylindrical topology). In this limit, the solutions of the field equations represent multi- 
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Skyrmionic configuration living in a flat tube-shaped region whose sections have a radius much bigger 
than the scale R. 

When N is large the equations read: 


0 = ^ -2(/ + l) j;(l + z)F, + 2iVj] +i^(l-cos(F0)F; 

V 1=0 k=0 i=k / 

- (1 - COS {Fi)) sin {Fi) 

+ [I (1 - cos (F,_0) + (/ + 2) (1 - cos (F,+i))] , (31) 

The first two terms in the above equation (involving summations) describe the free part of the theory, 
while the rest (involving R) is responsible for interactions, therefore defines nontrivial dynamics. 
Therefore we find our prescription (I30p a reasonable way to implement flat (large volume) limit within 
our formalism. In the next section we will demonstrate (for N = 3), that the dynamics does depend 
on R. 


5 Phase space portrait and geometrical ’t Hooft parameter 


One of the strongest advantages of our formalism lies in the fact that we have an autonomous system 
of differential equations. The latter can be qualitatively described by its phase portrait, an important 
characteristic of which is the set of critical points. In particular if under some change of parameters 
the set of critical points remains the same (their number and their character), then one does not 
expect that dynamics of the system changes dramatically under the same change. On the contrary, if 
the number and/or character of the critical points changes, one may definitely say that the dynamics 
of the system is qualitatively different. 

In this section we are going to study our system by using the tools of dynamical system theory 
(we will follow [431 144j 1. In particular, we will analyze how changing the value of R^ the number of 
the critical points of the dynamical system, as well as their linear stability properties, change. 

For simplicity, we shall perform the analysis in the SU (3) case. In this case we have two profiles, 
Fq and Fi, and the equations of motion are given by: 


Tn 


Fi 


2 I 

3 

2 I 

3 


(1 - cosFo) 


, I pi , sinTo 


(1 - cosFi) ) + + 5^ 


'^2 2(l-cosFb) , (1-cosTi) 

-4--+ 

f/ _ 4 _ + 




(1 - cosFq) 
i?2 


= 0 , 

= 0 . (32) 


Notice that the two equations go one into the other by exchanging Fq and Fi. We are going to exploit 
this symmetry. 

In order to study critical points of the system (1321) . hrst of all we rewrite it as a system of first 
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order equations as follows: 


where 


Z = 




A 



C(Eo,Fi,A) 


F 

, g(z) = 

Uj 

_ Q{Fi,Fo,E) _ 


Fo 

\ J 


A 


_ El _ 



E 



(33) 


(34) 


and we introduced the following notations for brevity 


a = 


2 1 1—cos Fq 

3 2 W 3 

1 2 I 1 1—cosFi 

3 3 2 W 




sin,^ 


r - 4- 


2(1 — cos^) (1 — cos??) 




+ 




(35) 

Stationary points Z* are defined by vanishing of the vector held G, which dehnes the how of the 
autonomous equation (133)) in phase space: 


(36) 


G (^Z* j = 0. 

One can easily check that deta > 0, therefore Z* can be written as follows: 


Z* — 


0 

0 

Fo* 

F* 


(37) 


where Fq and Fj^ are solntions of the following algebraic equations 

f Q(Fo*,Ff,0)=0 
1 Q(F*,F*,0)=0 ■ 


(38) 


Equations (I38p imply one of the three possibilities (the fourth possibility does not have real solu¬ 
tions): 

1. sinFo = 0; sinFi = 0; 

2. sinFi = 0; -4 - = 0; 


3. sinFo = 0; —4 — 


2(1—cosFi) (1—cos_Fo) _ 


'W 


~W 


= 0 . 


Linearizing (j33p in the vicinity of the critical point Z* we obtain 

- - = JE, € — Z - Z^:, 

ax 


(39) 
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where J is the Jacobian matrix of the diffeomorphism Z i—)• G{Z) taken at the point Z* 

, j,k = l,.A. (40) 

i=z* 

In order to analyze stability of the critical point Z* i.e. in order to understand whether phase trajec¬ 
tories starting near this point “run away” exponentially or not, one has to check whether eigenvalues 
of J have real part or they are purely imaginary. 

Below we present a complete analysis of the critical points. Since in the algebraic equations (j38p 
the unknowns Fq and Fi enter via sine and cosine, it is sufficient to consider critical points modulo 
2tt. 


T 

dZk 


5.1 Case 1 

We identify the following four independent subcases: 

.Fo* = 0,Ff = 0; 

• -^0 = 

.Fo* = 0,Ff = 7r; 

.Fo* = 7r,Ff = 0. 


5.1.1 Case = Ff = 0 

In the (equal prohle) case Fq = F* = 0 Jacobian matrix is 


J = 


0 0 
0 0 
I 0 
0 1 


2 

w 

1 

~w 

0 

0 


1 

~w 

2 

w 

0 

0 


(41) 


whose four eigenvalues are: 


Ai,2 = ± 


F’ 


As,4 = ± 


R ■ 


(42) 


Therefore this hxed point is unstable; this property does not depend on R. 
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5.1.2 Case = n 

Also in this case we have equal profiles Fq = F^ = tt. In this case the Jacobian matrix is 



0 0 

4 (2ij2+3j(^2ij2+ij 

1 2R?+1 



2 ~W(W+FWFW 

2 W+AW+3 



0 0 

1 2R?+1 

1 4iJ'‘+8_R2+S 


J = 

2 W+AW+3 

2 W{W+AW+T) 



1 0 

0 

0 



. 0 1 

0 

0 


and the eigenvalues 




^ V(«^ + 3)(2F^ + 1) ^ 

2 (F2 + 3)F ’ 

W2 v^(F2 + 1)(2F2 + 1) 
“ 2 (F2 + 1)F 


are purely imaginary, therefore this point is stable for all R. 


(43) 


(44) 


5.1.3 Case Fg* = 0, F^ = tt 

Now we consider Fq = 0, Fj* = tt. In this case the Jacobian matrix 

1 {2R^-l){2R^+3) 


J = 


0 0 
0 0 
I 0 
0 1 


R^{R^+2) 

1 2/?2-1 
2 W(W+T} 

0 

0 


WITF+T) 

2R‘^+l 

W(W+T} 

0 

0 


has the eigenvalues 


Ai,2 — ± 


As,4 — 4= 


(F2 + 2) VV48F8 + 120F6 + 24^4 - 48^2 + 9^^ 

ZR2 ’ 

i {R^ + 2) \/V48F8 + I20F6 + 24^^ - 48^2 + oTs 


(45) 


(46) 


Unlike the previous cases, now the properties of this critical point change at a critical value of the 
parameter R. Indeed, while the eigenvalues As and A 4 are always purely imaginary, Ai and A 2 are 
purely imaginary only when R < In this regime this critical point is stable. For R > I/V 2 

instead the first two eigenvalues become real, therefore the fixed point becomes unstable. This is a 
nice example of how the stability properties change at a critical value of the parameter. Interestingly, 
the value of R at which the transition happens is the same below which new fixed points appear (i.e. 
those occurring in cases 2 and 3 described above). 
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5.1.4 Case = vr, Ff = 0 


In this case the equations are exactly the same of case 3, only with 5 and e exchanged. The eigenvalues 
are therefore exactly the same so that the same discussion applies. 

5.2 Cases 2 and 3 

We now consider cases 2 and 3 in the above list. Since they are related by the exchange Fq ■(-)• Fi, they 
will have the same eigenvalues and hence the same stability properties. Let us then consider case 2: 


smn = 0 ; _4_2(l-™ro)^(l-cosn) ^ 


(47) 


F 2 F 2 

One may immediately check that a real solution only exists for F < 1/ y/2, and it is given by 

F* = TT, Fq* = arccos(2F2), (48) 

modulo 27r. For R > 11^/2 this critical point disappears. The Jacobian matrix reads 


J = 


44 44 (8^8-14^4+3) (2^2+3) 

^ ^ iJ2(4i?«-12iJ4+9) 

8fi8-14R4+3 

4 i?»- 12 R 4_|_9 


0 0 
1 0 
0 1 


_i 6i;2+3 

2 2R4-3 

1 12ij4-12fl;2_9 

'4 _R^(2i?4_3) 

0 

0 


(49) 


and its eigenvalues are: 

^ Y^(2F4 - 3) (2F2 + 1) (16F4 + 10F2 - 3 + + 896- 28AR^ - 924^2 + 441) 


Ai,2 = ± 


-^3,4 = ± 


(50) 


4 (2 F4 _ 3) 

^ yj- (2 F4 - 3) (2 F2 + 1) (-16 -10R‘^ + ^256 + 896 - 284 - 924 + 44T + 3) 

I (2F4-3)F 


In the whole range F < l/\/2, the first two eigenvalues are purely imaginary, while A 3 and A 4 are 
real, therefore this critical point is unstable. At F = l/\/2, A 3 and A 4 vanish. 

5.3 Summary 

The analysis we just performed discloses the crucial dependence of the dynamics of our system on 
the parameter F. In fact we saw that there is a critical value F = I/\/2 at which the nature and 
the number itself of the critical points change. For 0 < F < l/\/2 and for F > l/\/2 we have two 
qualitatively different phase portraits, thereby different dynamics. These are genuine finite volume 
effects (since the parameter F represents the available surface per Baryon within the tube-shaped 
region) which can be disclosed only within the present framework. 

In FiglUwe present location of critical points for F < I/\/2 and F > I/\/2. As one can easily 
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see, the two plots are qualitatively different: while the latter grid is invariant under the discrete group 
of reflections x about the axes Fq and Fi, the former is not, that is the symmetry is broken. 
Moreover the former is sensitive to the particular value of R, while the latter is always the same for all 
R > \l\pl. Finally, the concentration of different critical points changes as well, since some of them 
changed color (i.e. type) when R passed through the critical value, while new (green) fixed points 
appearecQ. 

These metamorphoses of the phase portrait hint at the presence of some sort of phase transition 
occurring at strong coupling (recall that XjR? plays the role of a coupling constant). It is interesting 
to notice that, in our units, this critical value is of the order of the /m, so although irrelevant from a 
flat limit perspective, it may have important consequences in view of applications to nuclear physics, 
where the dimensions of the tube-shaped region should be precisely of that order of magnitude in 
order to model an atomic nucleus |30] . 




Figure 1: Location of critical points for R < l/\/2 on the left and for R > on the right. In 

the latter case the grid is independent on a particular choice of R, while for the former case we chose 
R = Blue points represent stable fixed points, red points represent “absolutely” unstable fixed 

points (i.e. all eigenvalues of the corresponding Jacobi matrix are real), while green points represent 
“softer” unstable fixed points (i.e two eigenvalues values are real and two imaginary). 


^For the color code see the caption. 
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6 Numerical Solutions: the SU{4:) case 


In this section we perform a numerical study of the system in the = 4 case. In this case we have 
three profiles, Fq, Fi,F 2 , and the equations of motion reduce to: 


.^0(1 + (1 ~ cosFq) I + 


Rl 


2Fi + F 2 , sinFo 


+ 


Rl 


r ■ 2 

Fn -4- 


6(1 — cosFq) 4(1 — cosFi) 


Rl 


+ 


Rl 


= 0 , 


f ^ ^ - cosFi)^ + 


+ 


sinFi 

w 


^ 8(1 — cosFi) 3(1 — cosFq) 3(1 — COSF 2 ) 

1 4 : ^^75 F F 


Rl 


^2 ( 1 + ^ (1 - COSF2) ) F 

Rq 


Rl 

Fq F 2 Fi sinF2 
3 


Rl 


= 0 , 


r „ 2 6 (1 _ C 0 SF 2 ) 4(1- cosFi) 

-C 2 — 4- 7:^5 -h 


Rl 


Rl 


= 0 . 


The hrst and the second go one onto the other upon exchanging Fq o F2. We can exploit this 
symmetry to look for solutions which have Fq = F2 = F (another possibility allowed by the symmetry 
is to set Fq = —F 2 and at the same time Fi = 0, however this would give vanishing winding number, 
so we do not consider it). Let us also call Fi = G. In this case we get the two independent equations: 


2 \ F F 2G sinF 

^ ( I + ;^ (I - cosF) j F-^-F 


Rl 


^2 6 (l —cosF) 4(l —cosG) 

F - 4-^- - F - 


Rl 


Rl 


G ( I F (I - cosG)^ F F F 


RA 


^2 ^ 8 (I — cosG) ^ 6 (I — cosF) 


Rl 


Rl 


= 0, 
= 0, 


and the energy density (which is defined by F = FQ( 67 r) ^ f Sdx) is given by: 


.2 .2 1 

£ = 2F fG F2FGF-^ 

Rf 


•0 L 


3 F F 4 (1 - cosF) F 2 G F 4 (1 - cosG) 


F 


1 r 


F 9(1 — cosF) — 12 (1 — cosF) (1 — cosG) F 8 (1 — cosG)^ 

Rq 


(51) 


Before solving the system we must choose a set of boundary conditions to impose on F and G in such a 
way to get solutions with a nontrivial winding number. Recall that we have to impose (anti-)periodic 
boundary conditions on the held U: 


U 




(52) 
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For simplicity we shall limit ourselves to the periodic case. According to the ansatz we are using the 
field can be written as: 

U = exp jiao (^Pq “ ^ ~ ~ ^ 

= + AoPo) + AiPi) + A 2 P 2 ) • (53) 

It is convenient to use the second form. We notice that the matrix U is diagonal in the basis 
{|Fo)) 1 ^ 2 );}) therefore the above boundary condition is equivalent to its diagonal matrix ele¬ 

ments. The diagonal matrix elements of U are 

(Fo|C/|V"o) =et(2F+G)^ ^V,\U\Vi) = e-2^ = i{V2\U\V2)r. (54) 

Thus the boundary conditions are satisfied if 

F = F +2m7r, m e Z, G = G + dnvr, n G Z. (55) 

With these conditions the winding number is given by 

W = 6m + 8n. (56) 

We also need to impose conditions on the first derivative of U, i.e. of F and G. We choose the simplest 
possibility which is compatible with the periodicity of the first derivative of U, i.e. 



We have solved numerically the system (j5ip with the boundary conditions (j55p and (j57p . for some 
values of the integers m and n, taking Rq = L = 1, and we computed the energies of the solutions. 
The results are summarized in the tables below and in the one in sect. EH We have allowed m and 
n to vary in a small range, beyond which our numerical procedure is not very accurate. 



m 

n 

W 

E 

m 

n 

w 

E 

m 

n 

w 

E 

-1 

1 

2 

26.68 

1 

0 

6 

11.43 

1 

1 

14 

43.87 

-2 

2 

4 

100.65 

2 

0 

12 

42.41 

2 

2 

28 

168.31 

-3 

3 

6 

228.42 

3 

0 

18 

92.08 

3 

3 

42 

376.07 





4 

0 

24 

161.61 






The energies are exactly the same if we change sign simultaneously to m and n. There are sectors 
in which, to a very good approximatioij^l, the energy grows like the square of the modulus of the 
winding number, with different coefficients in different sectors. However, this is not true in all sectors, 
for example: 

®the accuracy is a bit lower in the last case of the third column, most likely because of numerical errors. 
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Figure 3: F,G and £ for m = —2, n = 2. 

m n W E 

0 1 8 24.48 

0 2 16 93.20 

We notice that the solutions with the smallest winding numbers are not energetically favored, since 
to achieve small winding numbers both profiles must grow going along x in opposite directions, while 
the smallest energy solutions are those where only one of the profile grows while the other oscillates 
and stays small, i.e. either m or n are different from zero, but not both. Notice also that increasing 
n costs much more than increasing m, despite their contributions to the winding number being not so 
different. 

In Figs.1-6 the plots of some solutions are reported. 





X 


Figure 4: F,G and £ for m = 1, n = 0. 
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Figure 5: F,G and £ for m = 2, n = 0. 



Figure 6: F,G and £ for m = 1, n = 1. 
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6.1 Comparison with the equal profile situation 


If we impose that the profiles are equal, i.e F = G = F, the equations of motion reduce of course to 
(j27l) . In this case the same procedure followed above tells us that the boundary conditions to achieve 
a periodic U are 

F ^ m' G Z (58) 

in which case the winding number is simply 


^^equal — 20 TTl . 


(59) 


In order to compare the equal profile situation with the non equal profile one we have to consider 
cases with the same winding number. This is achieved if in the latter case we have m = 2re. We did 
the computation in the two simplest cases, with Rq = L = 1, and the results are: 


m' 

equal 

Eequal 

m 

n 

W 

E 

1 

20 

82.67 

2 

1 

20 

82.67 

2 

40 

325.70 

4 

2 

40 

325.57 


We see that the energies are the same (in the second case the discrepancy is due to small numerical 
inaccuracies: the numerical computation for the m = 4, n = 2 case gave the winding number of the 
solutions to be 40.13, thus signalling such inaccuracies). Therefore, in topological sectors whose wind¬ 
ing number is compatible with equal profile configurations, with the chosen values of the parameters, 
equal-prohle solutions and non-equal profile ones are energetically equivalent. 


6.2 Dependence on Rq 

Besides the above analysis, which was performed at fixed Rq, we also studied what happens when this 
parameter is varied, considering for definiteness the lowest winding number case, i.e m = — l,n = 1. 
We checked that when Rq is increased the solutions tend to linear functions, in agreement with 
triviality of the dynamics for large radius. Instead, for values of Rq (or Rq/L if we allow also L to 
vary) substantially smaller than 1 the numerical procedure becomes unstable and the results unreliable. 
This is of course to be expected since lowering Rq the non-linearity of the problem increases. This is 
in agreement with the bounds (I89II94I) . which are more difficult to satisfy when Rq is small, and trivial 
when Rq is large. 

In the small Rq case a more sophisticated numerical analysis would therefore be needed. In 
particular it would be interesting to see whether at strong coupling non-equal profile solutions are 
energetically favored with respect to equal profile ones with the same winding number, as suggested 
by the fixed point analysis of section [5j Another interesting issue requiring a more accurate numerical 
work is the identification of a range of the parameters in which the dependence of the energy on the 
winding number is linear rather than quadratic. We hope to come back on these interesting issues in 
a future investigation. 
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7 Conclusions 


In this paper the SU(N) Skyrme model is studied in a tube-shaped geometry which allows to consider 
finite volume situations. The use of such a geometry allows to introduce a parameter which regulates 
the dynamics as a new (inverse) coupling constant. These simplihcations, which make the system 
autonomous, allow us to perform a detailed study of the system using several techniques, i.e. dynamical 
systems, numerical solutions and rigorous existence theorems. All these analyses confirm the relevance 
of the parameter Rq introduced by the metric on the dynamics of the system. As a further bonus, when 
combined with N as 1/R^ = N/{Rq)'^, this parameter defines an effective ’t Hooft coupling. In the 
large N limit in which such coupling is kept fixed, flat space-time is recovered. This effective’t Hooft 
parameter determines the discrete symmetries of the Skyrmions configurations. These results clearly 
show that the present framework is a very promising tool to study interacting Skyrmions systems. A 
very interesting problem which we hope to analyse in a future publication is to extend the present 
analysis at finite temperature and/or chemical potential. 
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Appendix: Existence of Solutions 

In this appendix we prove analytically that non-trivial solutions which are not trivial embeddings 
of SU{2) into SU{N) (in which the prohles are not proportional) do indeed exist. We will focus 
for simplicity on the SU{3) case but the same argument can be easily extended to the general case. 
The basic mathematical tool of this section is a well-known result of nonlinear functional analysis, the 
Schauder theorem (see for a detailed pedagogical review on mathematical tools to deal with non-linear 
partial differential equations |45] ). 

The statement of the Schauder theorem ( [351SS] ) is the following. Let S' be a Banach spac Let 
C a bounded closed convex set in S, and let T be a compact operatoi0 from the Banach space S into 
itself such that T maps C into itself: 

T[.]:C^C . (60) 

® A Banach space is a linear space endowed with a norm, and which is complete with respect to the metric induced by 
the norm. Recall that a metric space is a space in which a distance d {X, Y) between any pair of elements of the space 
is defined, and it is called complete if (with respect to the chosen metric) from every Cauchy sequence one can extract 
a convergent subsequence (see, for instance, [35]). 

^An operator T from a Banach space S into itself (see, for a detailed discussion, [451146] 1 is called compact if and only 
if, for any bounded sequence {xn}, the sequence {T{xn)} has a convergent subsequence. 
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Then the map T [.] has (at least) one fixed point in C. In other words, under the above hypotheses, 
there exists a solution to the equation 

T[X]=X . (61) 

Let us recall the equations of motion for the SU (3) case: 




1 


1 


^^'3^ 2R^ 


(1 - cosFo) 
(1 — cosFi) 


Ip sinFo 
+ 3^' + 

1 ” sinFi 
+ 3^° + 


2(l-cosFo) (1-cosFi) 

Fa -4-jp-+ 

F," - 4 - + 


F2 


F2 

(1 - cosFq) 
F2 


= 0, 
= 0 . 


(62) 

(63) 


We want to prove that the system in Eqs. (1621) and (16311 admits non-trivial solutions in which the 
profiles Fq and Fi are not proportional. In order to achieve this goal, let us rewrite it as coupled 
integral equations (for notational simplicity in this section we shall consider x € [0, L] instead of 
X G [—L/2,L/2], since everything just depends on the length L of the tube): 


where 


px 

Fq {x) = QQ + bQX- / Hq (s) ds , 
Jo 

px 

Fi (x) = oi + bix — Hi (s) ds , 

Jo 


Hois) 

Hi is) 

ho ip) 
hi ip) 
Ai 


f 

f 

Jo 


ho ip) + ip) > dp 


f sinFo ip) 

4ii I —(p)+‘'p 


'^2 2 (1 - cosFo (p)) , (l-cosFi(p)) 

Fo -4-jp-+-- 

^ ^ ^ 2(1- cosFi jp)) (1 - cosFo ip)) 


F2 


F2 


-1 


i = 0,1 


(64) 

(65) 

( 66 ) 

(67) 

( 68 ) 

(69) 

(70) 


where Fj = ■^, i = 0, 1 , and ai and bi represent the initial data for the two profiles Fq and Fi and their 
derivatives at x = 0. It is a trivial computation to show that the above system of integral equation is 
equivalent to the system in Eqs. (1621) and (1631) . The system in Eqs. (IMI) and (l6^ can be written as 
a fixed point condition for the following vectorial operator ^ acting component-wise on pairs of 
functions ^ (x) = (Fq (x) , Fi (x)): 

f : [0, L] X [0, L] [0, L] x [0, L] , 

f (x) = (Fo (x), Fi (x)) G [0, L] X [0, L] , 


^ [Fo,Fi] = 7 ^ [ (x)l = (To (x) ,Ti (x)) 


(71) 
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with 


To{x) 

Ti{x) 


ao + box 


ai + 61a; 



Hi (s) ds 


(72) 

(73) 


where the functions Hi (s) are defined in Eqs. (166)) . (1671) . (f68]l and (f69]l . It is then obvious that the 
system in Eqs. (|64p and (1651) can be written as the following fixed-point condition 


'f {x) = ^\^ 


X 


(74) 


where the operator T has been defined in Eqs. (|7T]) . (1721) and ([73]) . Hence, now the task is to define 
^ as a compact operator from a bounded closed convex sub-set of a Banach space into itself. 

In order to achieve this goal, first of all let us define the following metric into the space [0, L] x 

[0,L]: 


X)\+ sup 

x€[0,L] 


d(F,G)= sup \Fo{x)-Go{ 

^ ^ xG[0,L] 

-I- sup |Fi(x) - Gi(x)|-I- sup 

xG\^.L\ xG\0.L 


dFo 

dGo 

+ sup 
xe[o,L] 

d^Fo 

d?Go 

dx 

dx 

dx"^ 

dx^ 

dFi 

dGi 

-1- sup 

xG[0,L] 

d?F^ 

d^Gi 

dx 

dx 

dx^ 

dx^ 


xG[0,L] 3::G[0,L] dx xG[ 0,L] dx‘^ dx^^ 

f {x) = {Fo {x ), El (x)) ,t{x) = {Go (x ), Gi (x)) G [0, L] x [0, L] . 


(75) 


With respect to this metric, which is induced by a norm, the space [0, L] x [0, L] is a Banach 
space which we will call S. 

The next task to apply the Schauder theorem is to define a bounded closed convex sub-set G of 
the Banach space defined above (using the metric in Eq. (1751) 1 such that ^ maps C into itself. Let 
us define C as 


{:?(x) = (Eo(x),Ei(x)) G E| VxG [0,L] |Eo(x) - ao| < B , |Ei(x) - ai| < H , (76) 


1^1 

dx 


<B, 


d^Fo 


dx"^ 


<B, 


dFi 


dx 


<B , 


d^Fi 


dx"^ 


< B 


B G 


where Oj are the initial data appearing in Eqs. (I64p and ()65p so that G is closed by definition. It is 
easy to see that C is bounded as well, since from the definition of G it follows that 


|Eo(x)| <B + |ao| 
|Ei(x)| <B + |ai| 


dFo 

< B , 

d^Fo 

dx 


dx^ 

dFi 

<B , 

d^Ei 

dx 

dx^ 


<B , 


<B , 


(77) 

(78) 


which clearly implies that there exists M such that d{P < M, given any two functions ^ G C. 
It remains to prove that G is convex, i.e. to check that if (x) and ^ (x) both belong to G then also 
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ef (x) + (1 — 0) ^ (x) belongs to (7, V 0 G [0,1]. This can be seen as follows: 

\9Fi{x) + il-9)Gi (x)-ai| < \9 {Fi (x) - ai)\ + \{1 - 9) (Gi (x)-ai)| < 9B + {1 - 9) B < B , 

(79) 

\9Fo{x) + {l-9)Go (x)-aol < |0 (Fq (x) - ao)| + |(1 - 0) (Go (x) - ao)\ < 9B + {1 - 9) B < B , 

(80) 


while the conditions on the derivatives are trivially satisfied. 

Thus, C is closed, bounded and convex. The requirement that ^ maps G into itself will impose 
some constraints on the parameters, as we shall see shortly. 

From the definition of G in Eq. (I76p we deduce the inequalities 


sup \F[q (x)| < L 

3;G[0,L] 

sup |i7i (x)l < L 

xG[0,L] 


B 3 
'2 ^ ^ 

B 3 
'2 ^ ^ 


+ 4 + 


F2 


B + 4+^2 


(81) 

(82) 


sup 

x£[Q,L] 

sup 

a;G[0,L] 


dHo 

dx 

B 3 

- 2 8F2 


dHi 

dx 

B 3 

- 2 ^ 8F2 



(83) 

(84) 


(where we used the fact that Ai < 3/2), which will be needed later. 

The next task to apply the Schauder theorem is to prove that ^ is a compact operator. To do 
this we advocate the Ascoli-Arzela theorem, which states that if a sequence of functions (defined on a 
compact metric space) { F} is uniformly bounded and equicontinuous, then a convergent subsequence 
can be extracted from ilo. Thus compactness of ^ is equivalent to the statement that if '^n (x) = 


{Fq (x) , F” (x)) is a sequence in G, then the sequence ^ '^n {■ 

bounded and equicontinuous. 

Uniform boundedness can be proved as follows: 


X 


= (To (n; x), Ti (n; x)) is uniformly 


To (n;x)| 


ao + box 



< |ao| + |6oF| + F sup |Fo(n;x)| 
x€[0,L] 


and therefore, from Eq. (I8ip . one gets: 


To (n; x)| < |ao| + \boL\ + 


B 3 
'2 ^ ^ 


F2 + 4 + 


F2 


(85) 


^Recall that the sequence {fn} is uniformly bounded if \fn \ < M, where M does not depend on n; {/n} is said to be 
equicontinuous if, given e > 0, 3 5 > 0 such that |/„ (x) — fn (j/)| < e whenever \x — y\ <5 and, moreover, S does not 
depend on n (otherwise the sequence would be continuous but not equicontinuous: see [45] 1. 
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and, similarly, for |Ti (n;x)|: 


|Ti (n;x)| < |ai| + |6iL| + 


B 3 


B + 4+^2 


( 86 ) 


Now, the sequence of images ^ (x) has to belong to C as well. As it always happens (see 

and m) this will give some constraints on the range of the initial data as well as on the parameters 
B, L and R since, from the definition of C, one has to require that for any n 


\To (n;x) 
|Ti (n;x) 


ao\< B , 

dTo (n; x) 
dx 

= \bo - Ho (n;x)| < B , 

d^To (n; x) 
dx"^ 

= 

dHo (n; x) 
dx 

oi < B, 

dTi (n; x) 
dx 

= \bi - Hi {n;x)\ < B , 

d'^Ti (n; x) 

= 

dHi (n; x) 
dx 


<B, 
< B . 


(87) 

( 88 ) 


Consequently, as it can be easily seen comparing Eqs. (f87)l . ([88]) . (1761) . ([771) and ([78|) with Eqs. (ISD) . 
([82|) . (f83]l . (I8^ . (I85|) and (l86t) . the following constraints arise; 


|6oA| + 


B 


+ 


\b,L\ + L^{- + 


2 8i?2 

B 


2 8i?2 


5^ + 4 + 


5^ + 4 + 


A 

6 


<B, 


<B , 




B^ + i + 


5 ^ + 4 + 


6 

6 


<B , 


<B , 


(89) 

(90) 


(91) 

(92) 


B 3 

B 3 




<B , 
< B . 


(93) 

(94) 


Therefore, in order for this theorem to work, the length L of the tube-shaped region in which these 
multi-Skyrmions are living cannot exceed the bounds defined in Eqs. (I89p and (I90j) (indeed, if L is 
too large, the bounds will be violated at a certain poinl]^.). It is also to be noticed that one cannot 
obtain a very large value for the allowed L by increasing B since the left hand sides of Eqs. (|89p and 
(|90p increase faster than the right hand sides. Moreover, the situation gets worse if is very small 
(as all the above inequalities are violated at a certain point). However, if R is large (namely, in the 
flat limit), Eqs. ([93P and (IMl) are always satisfied and Eqs. (IMP and ([MP become mild constraints on 
the initial data. On the other hand, it is trivial to see that it is always possible to choose the initial 

®the fact that L cannot be arbitrarily big for the theorem to hold can be also seen by observing that in the limit 
L —^ 00 the domain of definition of the functions would not be compact any more, thus invalidating the Ascoli-Arzela 
theorem 


25 












































data and B, L and R in snch a way that all the above ineqnalities are fnlfilled. 

The second step to prove that T is compact is to show that if (x) = (F-(x),F-(x)) i 


IS a 


seqnence in C then the sequence ^ (x) is equicontinuous. To show this, we must evaluate, for 


a generic n, the absolute values of following differences: 


\To {n;x) - Tq (n;y)| = 


|Ti (n-x) - Ti {n;y)\ = 


ry 

bo{x-y)+ / Hq (re; s) ds 
J X 
ry 

bi{x - y) + / Hi (re; s) ds 
J X 


(95) 

(96) 


where 0 < x < y < L. After trivial manipulations (which use the fact that all the functions '^n (x) 
belong to C and consequently Eas. (f77D and (1751) are satisfied) one arrives at 


\To (re;x) - Tq (re;y)| < |x - y| 
|Ti (re;x) - Ti (re;y)| < |x - y| 
Thus, given any e > 0, one can choose 
h < 




B +4+^2 


5^ + 4 + 


B? 


‘^{\bo\ + \bi\ + L{^B + ^[B‘^ + 4+^]}} 
in such a way that both the choice of 5 in Eq. (|99l) does not depend on re and, moreover. 


(97) 

(98) 

(99) 


\x — y\ < 5 


In summary, Eqs. 


\Ti (re; x) -Ti{n;y)\ < - , V re, i = 0,1 . 

, ([TO . ([Mil . (l92l) . (fM|) and (IMl) show that, if 


( 100 ) 


X = 


is uniformly bounded in C. 


X 


in C, then the sequence ^ (, 

theorem, from any sequence ^ 


X 


{Fq (x) ,Ff (x)) is any sequence in C, then the sequence ^ n {x) 

Secondly, Eqs. (f?7l) . (fM]l . (l99)l and (llOOp show that, if '^n{x) = {Fq (x), Ff (x)) is any sequence 

is equicontinuous. Consequently, by virtue of the Ascoli-Arzela 

one can extract a convergent subsequence: this, together 
with the bounds (I89|) - (|94|) . implies that the operator ^ is a compact operator from a bounded closed 
convex set into itself. Thus, it is possible to apply the Schauder theorem, which ensures that Eq. (1711) 
(which is equivalent to our original system) has at least one solution. Moreover, it is always possible to 
choose appropriately the inital data a* and bi in such a way that the two profiles are not proportional. 

The conclusion is that not only one can construct numerical solutions (which are very interest¬ 
ing by themselves), as we do in the main text, but also one can prove analytically that non-trivial 
multi-Skyrmions in which the profiles are not proportional do indeed exist. Besides the intrinsic math¬ 
ematical elegance of the fixed-point Schauder-type argument, the present procedure also discloses the 
presence of the bounds in Eas. (|89p - (l94p on the radius R and on the length L of the tube-shaped region 
in which these multi-Skyrmions are living (which constrain its shape), as well as on the other param- 
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eters of the model. At the present stage of the analysis, it is not possible yet to say whether such a 
bound is just a limitation of the method or it signals some deeper physical limitation on the volume 
of the regions in which one constrains these Skyrmions to live. However, to understand whether or 
not SU {N) multi-Skyrmions can fit into very large tube-shaped regions is certainly a very interesting 
and deep question on which we hope to come back in a future investigation. 
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